Coulomb wave function correction including 
momentum resolution for charged hadron pairs: 

Analysis of data for tt + tt^ pair production in p + Ta reaction 

at 70 GeV/c 

T. Mizoguchi 1 : M. Biyajima 2 } I. V. Andreev 3 * and G. Wilk 4 § 

1 Toba National College of Maritime Technology, Toba 51 7, Japan 

2 Department of Physics, Faculty of Science, 
Shinshu University, Matsumoto 390, Japan 

3 P. N. Lebedev Institute of Physics (FIAN), 117 924 Moscow, Russia 

4 Soltan Institute for Nuclear Studies, Nuclear Theory Department 
Hoza 69, PL-00-681 Warsaw, Poland 

Preprint DPSU-96-9 and SINS-PVIII-1996-15 

February 1, 2008 



Abstract 

We propose a new method for the Coulomb wave function correction including 
the momentum resolution for charged hadron pairs and apply it to the precise data 
on ir + 7r~ correlations obtained in p + Ta reaction at 70 GeV/c. It is found that 
interaction regions of this reaction (assuming Gaussian source function) are 5.6±3.0 
and 4.4 ± 2.6 fm for the thicknesses of the target 8 and 1.4 microns, respectively. 
The physical picture of the source size obtained in this way is discussed. 
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1. Introduction. Recently we have obtained the new formulae for the Coulomb wave 
function correction for charged hadron pairs ]T|, fj]. In particular we have applied them 
(in [Q) to data on tt + tt~ correlation obtained in p + Ta reaction at 70 GeV/c || (which 
were originally corrected by usual Gamow factor only). However, as it was pointed out 
to us by one of the author of ||, we did not take into account their published finite mo- 
mentum resolution [|J. In fact, our formulae cannot be applied directly to experimental 
data in which such momentum resolution is accounted for. Therefore in the present letter 
we would like to extend our method for the Coulomb wave correction provided in |], 0] 
to include also the momentum resolution case and to re-analyse data of || and also to 
analyse the new, preliminary data of || obtained with two kinds of thickness of Ta target: 
8 microns (8 mkm) and 1.4 microns (1.4 mkm). 

In the next paragraph we first reconstruct (for the sake of completeness) the analysis 
performed in Ref. 0] and then, in the third paragraph, we propose our new method for 
the Coulomb wave correction including this time also the momentum resolution. The 
final part contains our concluding remarks. 

2. Reconstruction of the analysis of 7r + 7r correlation data performed by using 
the Gamow factor with momentum resolution. It has been stressed in || that 
relative momenta of 7r + 7r~ pairs observed by them have some finite resolutions. The 
averaged correlation function, defined as 

R( + ( * '^' I [ * ^° * ^° (1 + B ) \ \ (1) 

\adpidp 2 / \adp 1 adp 2 w ) / ' 

depends therefore on this momentum resolution, where (1 + B pr ) stands for non-Coulomb 
correlation factor. To account for it the following random number method has been pro- 
posed in Ref. || in order to obtain the corresponding averaged quantities in analysis of 
the correlation data. 

First of all, the relative momentum of the measured pair, q = pi — p 2 , is decomposed 
into its longitudinal and transverse components: q L and q T , respectively, by making use 
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of the uniform random number u G (0, 1) (it is worthwhile to notice at this point that 
the transverse components q T in data of || are smaller than 10 MeV/c). One uses the 
following scheme here: 

10y/u for q > lOMeV/c, , 
q T = { _ (2) 

q^fu for q < lOMeV/c , 



<l \/q 2 - Qt ■ 

At the next step, the Gaussian random numbers for q L and q T are generated in the 
following way: 

<?L(random) = CTl-^ + ^l, (3) 
9T(random) 

X + q T . (4) 

In the above equations X stands for the standard Gaussian random numberf] whereas cr L 
and cr T are longitudinal and transverse setup resolutions for the corresponding compo- 
nents, which are equal to (values used in 0): cr L = 1.3 MeV/c; a T = 0.6 MeV/c (for 
target of the thickness 8 mkm) and cr T = 0.4 MeV/c (for the 1.4 mkm target). 



Finally, using the randomized number g random = ^(random) + ^(random) one calculates 
the corresponding randomized Gamow factor correction: 

\ 27rr] ran ij om . . 

l — "^random ) — , „ \ 1 ' V / 

exp(-27r?7 random ) - 1 

where ty-andom = m «/<?random- The full flow chart for this procedure is shown in Fig. 1. 
Calculating now the average value of G{— r/random) in 100 k events one can estimate the 
Gamow factor with this finite momentum resolution, 

R( q ) = G(- V )b+(l-b), (6) 



1 The probability density of such random numbers is given by 

f(x) — j_ exp(— x 2 /2) (—oo < x < oo) 

V 27T 
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where b is a free parameter. It is understood (or, rather, implicitly assumed) that essen- 
tially all unlike sign pions one deals with here originate from decays of long lived particles 
like 77, Kq , A, and so on[|. Figure 2 shows the results of analysis of new data (for 8 mkm 
target) || for region q > 3 MeV/c using this method. 

3. Proposal of the new method. We would like to propose now a new method 
of Coulomb wave function correction with a source function p(r) instead of the Gamow 
factor, in order to analyse the same data. As usual we decompose the wave function of 
unlike charged bosons with momenta p\ and P2 into the wave function of the center-of- 
mass system (cm.) with total momentum P — \{pi + P2) an d the inner wave function 
with relative momentum q = {jp\ — p2)- This allows us to express Coulomb wave function 
^(q, r) in terms of the confluent hypergeo metric function $ j7|: 

#(q, r) = r(l - zr/)e^ /2 e iq - r/2 $(zr7; 1; iqr{\ - cos 0)/2) , (7) 

where r = x\ — Xi and the parameter rj = ma/q. Assuming factorization in the source 
functions, p(r 1; r 2 ) = p(ri)p(r 2 ) = p(R)p{r) (here R = \(x\ +x 2 )), we obtain the expres- 
sion for Coulomb correction for the system of 7r + 7r~ pairs identical (modulo the sign) as 
in IB!: 



C (-rj) = J p(R)d 3 Rj p(r)d 3 r\V(q,r) 



^ ^ n + m + l 

G (- V )[i + A lc (- V )} , 



where 



J R (n,m) =4tt J dr r 2+n+m p(r) , A n 



T(—ir] + n) 1 



r(— iff) (n\) 2 

For the specific choice of Gaussian source distribution, p(r) = -C= exp(— (3 2 r 2 ), we have 



2 fl\ n+m ^ fn + m + 3\ 



2 This fact has some consequences on the determination of the radius of the interaction region as will 
be shown later on and makes it different from that obtained from the Bose-Einstein correlations. 
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whereas exponential source function, p(r) = |^-exp(— /3r), leads to 

Using now the same method of Gaussian random numbers as in the previous paragraph, 
we can analyse the old and the new data on 7r + 7r~ pairs H, |5j using the following formula: 

R(q) = C c (-rj)b + (1 - b) . (11) 

Figs. 3 and 4 show results of our analysis of the old and the new data, respectively. Ta- 
ble | show our results obtained using eq. (H) applied to old and new data with q > 3MeV/c. 

4. Concluding remarks. We have proposed the new method for the Coulomb wave 
function correction with momentum resolution and applied it to the analysis of the precise 
data provided by [Q, ||. Authors of Ref. || have analysed their 7r + 7r~ correlation data using 
Gamow factor for Coulomb corrections together with the random numbers method to 
account for final momentum resolution. We have repeated this analysis replacing Gamow 
factor by the Coulomb wave function but following the same method for correction for 
the momentum resolution effect (cf. eq. (§)). As a result we have obtained the following 
ranges of interaction for the Gaussian source function: 

r(p + Ta) = — = 5.6 ± 3.0 fm for 8 mkm, (12) 
Zp 

= 4.4 ± 2.6 fm for 1.4 mkm. 

To get a correct physical picture of the source size, we should calculate the root mean 
squared size, which is equal to: 

a/3 

r rms = — — = 9.7 ± 5.3 fm for 8 mkm, (13) 
zp 

= 7.7±4.5fm for 1.4 mkm. 

The present study of tt + tt~ pair correlations has shown therefore that one can estimate 
the interaction region even from the 7r + 7r~ correlation data. It can be compared with the 
size of the Ta nucleus, which is given by: 

r Ta = 1.2 x A 1/3 

= 1.2 x (181) 1/3 = 6.8 fm. (14) 



As one can see, r rms is significantly bigger than r Ta . We attribute this difference to a 
physical picture shown in Fig. 5, i.e., to the fact that unlike-sign pions are mostly (if not 
totally) emerging from the long-lived resonances shown there. (In a future one should 
consider also a possibility of more direct estimation of the parameter b and its role in 
determining the source size parameter^). 

For completeness we have also tried to analyse the same data using exponential source 
function instead of Gaussian. As is shown in Table |I| this leads to errors on r = of the 
order of 100%, i.e., with this type of source function we cannot estimate the source size 
(therefore it has to be discarded). 
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3 In an analysis of 7r°7r° correlation data, a similar function f(q) is introduced: 

R(n\ a ) = f(q) + (l-f(q))ll + \E 2 2B }, 

where A and E 2 b are the degree of coherence and an exchange function due to the Bose-Einstein effect. 
f(q) is attributed the resonances effect; f(q) » 0.9 ~ 0.7 depends on the Monte Carlo programs 0. 
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Table I: Results of the \ 2 fits of R(q) for Gaussian source by eqs. (§) and (pTTf) . 



Reaction 


Formula 


1/2/5 [fm] 


b 


X 2 /NDF 


data of Ref. § 


Gamow factor 


— 


— 


57.8/40 


cf. Ref. @ 


eq. @) 


2.30 ±0.88 


— 


51.0/39 


(without momentum resolution) 










data of Ref. 


eq. © 


— 


b = 1 (fixed) 


53.9/37 




eq. org) 


2.96 ± 1.02 


b = 1 (fixed) 


44.9/36 


data of Ref. § 


eq. © 




0.43 ±0.03 


55.0/46 


8 mkm 


eq. © 


5.62 ±3.03 


0.53 ±0.07 


51.5/45 


data of Ref. § 


eq. © 




0.51 ±0.04 


37.9/46 


1.4 mkm 


eq. © 


4.44 ± 2.63 


0.60 ±0.07 


35.1/45 



Table II: Results of the \ 2 fits of R(q) for exponential source by eq. (|TT|). 





gum [MeV] 


1/(3 [fm] 


b 


X 2 /NDF 


8 mkm 


25 


3.19 ±3.55 


0.51 ±0.08 


24.0/20 


1.4 mkm 


25 


3.00 ±3.88 


0.59 ±0.10 


20.3/20 




30 


0.65 ±2.62 


0.54 ±0.08 


23.1/25 




35 


2.08 ±2.77 


0.57 ±0.08 


27.2/30 




40 


2.60 ±2.85 


0.58 ±0.08 


30.3/35 
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Figure Captions 

Fig. 1. Flow chart of the present procedure for Gamow factor and/or Coulomb wave 
function, which includes momentum resolution by generating the Gaussian random 
numbers. 

Fig. 2. Results of x 2 fit f° r 8 mkm target of p + Ta — > 7r + 7r~ + X reaction with q > 3 
MeV/c by eq. (§). 

Fig. 3. Results of the \ 2 fit for Ref. of p + Ta — > -n + it~ + X reaction with q > 3 
MeV/c by eq. flll]). 

Fig. 4. Results of the \ 2 fits f° r P + Ta — > 7r + 7r~ + X reactions with q > 3 MeV/c by 
eq. (|TT|) : (a) 8 mkm target; (b) 1.4 mkm target. 

Fig. 5. Physical picture of the obtained source size based on the assumed resonance 
effects. 
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Generation of uniform random 
numbers: u v u 2 



S = (2u 1 -l) 2 + (2u 2 -l) 2 




No 



X 1 = (2u 1 - 



X 2 = (2u 2 - 1) 



-21ns 



-21nS 



^L(random) 




^T(random) — 


: (5jX 2 + 



^random ^ ^L(random) ^T(random) 



Gamow factor / Coulomb wave function 



10 k ~ 100 k events 



Fig. 1 
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Fig. 2 
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Fig. 3 
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(a) Target 8 mkm of Ta 
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Fig. 4 
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(b) Target 1 .4 mkm of Ta 
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Fig. 4 
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R~1.2A 1/j ~6.8fm 



R„ m -11.8 fin 



Fig. 5 
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